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The heavy fermion compound URu2Si2 exhibits 
some of the most enigmatic manifestations of 
strongly correlated electronic matter. Cooling 
down to To — 17. 5K, reorders its ground state via 
a second order phase transition to an unidentified 
exotic order, the so called 'hidden order' (HO), 
whose fingerprint is a large specific heat jump and 
a disproportionally small antiferromagnetic mo- 
ment. Below 1.5K, the HO becomes supercon- 
ducting, while it evolves to a large moment anti- 
ferromagnetic state (LMAF) upon raising pres- 
sure, through a first-order phase transition. A 
cascade of metamagnetic steps, an anomalous gi- 
ant Nernst signal and a characteristic gapped 
spin excitation at Q = (1,0,0), add more mys- 
tery to the HO's identity. Here we provide for 
the first time, a unified understanding of the 
observed anomalous thermoelectricity, the unre- 
solved metamagnetic behaviour and the existence 
of two distinct energy scales, demonstrating that 
the HO state is in fact a chiral d^y -\- id^2_^2 spin 
density wave state. The inelastic magnetic re- 
sponse at Q, is attributed to the emergence of a 
novel orbital antiferromagnetic collective mode, 
rotating the HO to the LMAF. Finally, we make 
specific predictions for a Polar Kerr experiment. 

The 'hidden' order (HO) transition of URu2Si2 at- 
tracted enormous attention, due to the mismatch of its 
specific heat jump at To=17.5K and the concomitant c- 
axis oriented tiny antiferromagnetic moment 0.03/^s/[/, 
of wave-vector Qq = (0, 0, 1)^. This incompatibil- 
ity lias raised great controversy on wliether the genuine 
spin-dipolar antiferromagnetic moment inferred from the 
Bragg reflection at Q(j^, may be the driving order pa- 
rameter of this transition. In fact, NMR measurements 
have revealed that antiferromagnetism is inhomogeneous 
in the HO state and the small observed moment arises 
probably due to a reduced volume effect^i^. In addi- 
tion, experiments have found no evidence for higher spin- 
correlators^. These observations are strongly in favour of 
a multipolar order parameter scenario. 

Inelastic neutron scattering (INS) measurementa^rli, 
have revealed two gapped spin excitations at the wave- 
vectors Q = (1,0,0) (equivalent to Qq) and Qi = (1 ± 
0.4,0,0). Specifically, the first resonance always accom- 
panies the HO state and disappears in the LMAF— i^i, 
while the second persists and acquires a larger gap by 
effecting pressure. In the HO, the gap of the Q exci- 
tation, uIq, increases upon raising the field, while the 
ujq slightly diminishes^. The latter implies the exis- 



tence of two different energy scales, that none of them 
seems to fit to the energy gap dictating the thermody- 
namic properties of the HO^. Moreover, the emergence 
of incommensurate excitations can be only understood 
within an itinerant picture^ that is also supported by the 
existence of strong Fermi surface nesting at Qi— and 
the recent verification of Fermi surface gapping in An- 
gularly Resolved Photoemission Spectroscopy (ARPES) 
experiments^^. In fact, the coincidence of the nesting 
vector with the momentum of the incommensurate spin 
excitation has been consideredi^ as the consequence of 
the creation of a spin cxciton, due to the formation of a 
density wave state with ordering wave- vector Q^. 

All the aforementioned properties point to the direc- 
tion of a multipolar density wave state developed at 
the incommensurate vector Q^. The quadrupolar im- 
plementation of this scenario had been proposed ear- 
lier by Ramirez et afi^, in order to account for the 
anomaly observed in the non-linear susceptibility. Nev- 
ertheless, a proposal of this type has to answer first of 
all the emergence of tremendously important character- 
istic phenomena that accompany the HO phase. These 
arc the puzzling metamagnetisroiiii^ and giant Nernst 
signaliSi^S, together with the existence of two different 
energy scales^ii^. 

Here, we demonstrate that the hidden order that lurks 
in the B — T phase diagram of URu2Si2 is the incom- 
mensurate chiral dxy + id.j;2_y2 spin density wave state 

ordered at Q^. This state consists of a two component 
order parameter, with the second of them corresponding 
to the orbital antiferromagnetic orders proposed earlier 
for the HO^"— . The additional d^y component renders 
the HO chiral rather than orbital, violating time-reversal 
in a macroscopic level due to the development of an in- 
trinsic ferromagnetic orbital moment. It has been shown, 
that due to this orbital magnetic-field-coupling, any un- 
conventional density wave of the form idr^2_y2 will ac- 
quire an additional field-dependent d^y component^liS^, 
that could provide the missing energy scale. Even if the 
chiral order does not appear in zero-field, it will be neces- 
sarily field-induced affecting the B — T phase diagram of 
the HO and every transport measurement realized in the 
presence of a field, such as the Nernst effect. Notice that 
associating the HO with a chiral state may be crucial for 
understanding the development of an unconventional su- 
perconducting state confined in the HO, that has been 
also considered to be chiral^^. 

The proposed two-component chiral state may be 
viewed as the most general picture of an unconventional 
density wave order. In this sense, our theory is fully 
compatible with all the earlier proposed itinerant uncon- 
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ventional density wave models, and it could be consid- 
ered as a natural extension to them. However, the chi- 
ral state exhibits properties not familiar to orbital or- 
ders. The violation of time-reversal symmetry leads to 
a Polar Kerr effect^ and the generation of anomalous 
Hall electrio^^ and thermoelectric current a^^'^° , inclu- 
ding the 'tilted-hiir giant Nernst signal recently shown 
to arise in strongly insulating chiral condensates^. By 
adopting a simple microscopic model and an S0(2) in- 
variant Ginzburg-Landau theory, we shall explain in a 
unified manner not only the metamagnetic behaviour and 
the anomalous thermoelectric transport but also the ine- 
lastic neutron scattering spectra. 

Within a mean-field approximation, we have repro- 
duced sclf-consistcntly the B — T phase diagram of this 
compound. For mediating the HO we have considered 
repulsive nearest and next-nearest inter-site Coulomb in- 
teractions of the form "Ylii^j ^ij^i^j- "^^^ order parame- 
ter field is described as a c-axis polarized spin-triplet or- 
der parameter Aq^(A;) — Ai sinfc^ sinfcj, — i/S.2{coskx — 
cosfcj,), with ordering wave- vector Qi. For the metal- 
lic state of URu2Si2 we have considered a single-band 
nesting model described by the energy dispersion £{k) ~ 
— 2t(cos kx + cos ky). We have to remark that our results 
are independent of the microscopic details of the partic- 
ular model. 

The resulting B — T phase diagram is depicted in 
FiglT^,b. We observe that our phase diagram grasps to a 
large extent both qualitative and quantitative features of 
the experimental observation s . At S = OT, Ai = 
and A2 = 1.55meV. The latter order parameter corre- 
sponds to the idx2-y2 component that as we observe in 
FigHt, slightly decreases upon raising the magnetic field 
for low temperatures. In contrast, the dxy order param- 
eter Ai, rises monotonically. We locate a Metamagnetic 
Critical End Point^ (MCEP) at Bci = 33. 5T and T~3K, 
down to which the driving HO gap, A2, evolves contin- 
uously with the magnetic field satisfying the quadratic 
relation A2(Sz)/A2(0) ~ 1 - (BJB^^)^. This behaviour 
has been verified in magnetoresistance measurements^^. 
In stark contrast, A2 remains practically constant for 
temperatures below the MCEP, giving rise to the field 
independent energy scale witnessed in the inelastic mag- 
netic response at Qi—. A continuous evolution with the 
field is not compatible with the first-order character of 
the metamagnetic step. This property settles the contro- 
versy concerning the field dependence of the HO. 

Our self-consistent phase diagram exhibits double-stcp 
metamagnetism (FigHJl) at i3ci=33.5T and Bc2 c^41T, 
below 3K in agreement with experiment o^^'^^ . Notice 
that the occurrence of metamagnetism in such a high 
field scale, is intimately related to the appearance of the 
dxy component, leading to an enhancement of the criti- 
cal temperature of the dominant d^^^y^ order parameter. 
The gradual increment of the field, shifts the four quasi- 
particlc bands through the combined orbital-Zeeman cou- 
pling. At each Be-, an additional energy branch crosses 
the Fermi level (FiglT];), providing an abrupt increase in 
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FIG. 1: a. Self-consistently extracted B — T phase diagram of 
the chiral d-spin density wave. The chiral HO is present every- 
where in this phase diagram. The reduced magnitude of the 
order parameters in the high field and temperature regions, 
mask its presence providing a thermodynamic phase bound- 
ary (dashed lines), b. A Metamagnetic Critical End Point 
(MCEP) is located at T ~ 3K and B ^ Bd = 33. 5T. Two 
metamagnetic steps appear at Bci ~ 33. 5T and Bd — 41T 
along with the formation of crossover areas (separated by 
dash-dotted lines), in perfect agreement with the experimen- 
tal observationsi^ii^. For Bci < B < Bc2, the chiral HO may 
also coexist with other still unidentified field-induced states, 
c. For T<3K, the d^2_y2 order parameter is nearly constant 
until B = Bel, where it abruptly decreases. This is the gap 
seen in neutron scattering experiments^. Above the MCEP, 
the gap evolves smoothly with the field, explaining the de- 
pendence arising from bulk measurements^^. The saturation 
appearing in high fields does not give rise to any thermody- 
namic anomaly, rendering the chiral HO experimentally un- 
detectable, d. The metamagnetic steps are accompanied by 
an abrupt increase of the magnetization. For higher temper- 
atures the first order character becomes milder leading to a 
smooth evolution, e. The cascade of the metamagnetic tran- 
sitions, appears due to additional band crossings at each S^— . 
The reconstruction of the Fermi surface increases steeply the 
carriers of the system and concomitantly the magnetization. 



the density of carriers of the system and a concomitant 
steep enhancement of the magnetization^. For higher 
temperatures the first order character of the metamag- 
netic steps becomes milder, leading to a smoother in- 
crease of the magnetization in this region of fields. 

The pronounced similaritities presented in our phase 
diagram compared to Ref i^^i^^ , are accompanied by im- 
portant ramifications arising from the finite chirality of 
the state. First and more importantly, the chiral HO 
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is present everywhere in this diagram, due to the en- 
hanced critical temperature, originating from the orbital 
magnetic-field coupling. For magnetic fields Bci < B < 
Bc2, both order parameters of the chiral state suffer a 
severe decrease in their magnitude. Consequently, the 
chiral d-spin density wave state is still present in this re- 
gion but unavoidably masked due to the destruction of its 
topological rigidity. The HO persists even for B > Bc2, 
behaving as a polarized paramagnetic metal. 

Despite of the continuous presence of the chiral d-spin 
density wave order, the emergence of the MCEP leads to 
the generation of novel states at Bci, that in this case 
coexist with the chiral HO constituting phases ///, V in 
Refji^. It has been shown that, under quite generic con- 
ditions, the effect of a Zeeman field on a charge density 
wave can lead to a field-induced Spin Density Wave tran- 
sition and vise versa^SiS. In this sense, the magnetic field 
may permit a number of possible induced states^ such as 
the chiral spin singlet d-density wave state, studied in the 
context of the high-Tc superconductor a^^i^^'^^i'^°'^^ , or 
even the spin singlet and triplet Pomeranchuk states^i^. 

Apart from the phase diagram, we have also studied 
in detail the thermoelectric response of this state. Ac- 
cording to experimental observations, this compound is 
characterized by an anomalous thermoelectric behaviour 
up to 35 T^^i^° , where the HO has been considered to 
collapseSS. The arising Nernst signal has a positive sign 
and a giant magnitude of several fiV/K, while the tem- 
perature evolution of the Nernst response resembles to a 
'tilted-hiir, demonstrating a peak at about 3-4K, where 
it reaches the value of 30/iV/K for B ~ I2T. Up to now, 
the enhanced Nernst signal has been viewed as the signa- 
ture of the semi-metallic character of this electron-hole 
compensated meta l^^i^° . 

Here, by taking into account both topological and 
quasiparticle sources in computing the electric (cxx, o'xy) 
and thermoelectric conductivity {axx,C(xy) tensor ele- 
ments, we demonstrate that the Nernst response in the 
HO state can be understood in terms of the chiral d- 
spin density wave state. Using the values of the order 
parameters providing the phase diagram of FiglU we 
obtain a giant Nernst signal, that has a magnitude of 
about one order larger than the one observed in the ex- 
periments. The origin of the 'tilted-hill' Nernst signal 
is directly related to the chirality and the strong Fermi 
surface gapping. As presented in Refi^, the strong insu- 
lating character of the state forces the thermopower and 
the Nernst signal to become equal at a thermoelectric 
crossing point. In the vicinity of this point, the usually 
large values of the thermopower also imply a Nernst sig- 
nal of the same magnitude, leading unavoidably to an 
enhancement of the latter. However, the thermoelectric 
rigitity and the tendency towards a giant Nernst signal 
are lost exactly at Bci- The steep decrease of the d^^-ys 
gap across this first order transition and the concomitant 
Fermi surface reconstruction which partially destroys the 
insulating properties of the state, lead to the full supres- 
sion of the response (Figl2^) in accordance with Refi^. 
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FIG. 2: a.. B — T dependence of the self-consistently calcu- 
lated Nernst signal. The arising Nernst response presents a 
large positive signal, of about one order larger than the exper- 
imentally observed^. The giant response originates from the 
topological robustness of the chiral stated. The field destruc- 
tion of the Nernst signal at Bci agrees with the experimenl^^^ 
although here, it is accompanied by the dramatic reduction, 
rather than the collapse of the dxy,d^i_y2 order parameters, 
b. The Nernst coeflicient u — N/Bz shows a 'tilted-hill' pro- 
file, with a pronounced peak in the region 3 — 5K following the 
experimental datc»i^. c. Field dependence of the resistivity 
for several temperatures. For T < To, the presence of a 'kink', 
also found in experiments^^, may be naturally attributed to 
the field-enhancement of the dxy component. 



Furthermore, the Nernst coefficient exhibits a peaked 
structure with temperature and decreases upon raising 
the magnetic field already at lOK (Figl^b). Notice that 
although this state is also present even above 17. 5K, it 
does not lead to a large Nernst signal, becoming in this 
manner experimentally elusive. We have also calculated 
the resistivity for finite fields and temperatures. We ob- 
serve in Figl2t that the field dependence of the resistivi- 
ty exhibits a kink with the increase of the magnetic field. 
This is a definite fingerprint of the field induced dxy com- 
ponent and coincides with the experimental findings^S. 
This kink shifts to lower magnetic fields upon raising 
temperature, since the dxy gap gradually weakens. 

The HO exhibits an equally intricate behaviour un- 
der the application of pressure. Specifically, it switches 
to a large moment antiferromagnctic state (LMAF) of 
wave-vector Qq = (0,0, 1), for pressure Pc = O.SGPa at 
T=OK— lii. The intimate connection between the HO 
and LMAF— and the emergent adiabatic continuity^, 
has motivated the authors of Refi^i^ to propose a 
U(1)=S0(2) theory for the imification of the HO and 
LMAF. We demonstrate that within our model, invoking 
an SO (2) symmetry can fully explain the INS spectra. 

In our case the dx2-y2 order parameter is the dominant 
chiral HO component while the dxy part is field-induced 
and depends on the former. As a matter of fact, the 
SO (2) symmetry should emerge between the idx2_y2 spin 
density wave (dSDW) state of wave- vector Qi = (I ± 
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FIG. 3: a. The p — B phase diagram obtained by the S0(2) 
invariant Ginzburg-Landau model employed for unifying the 
hidden order (HO) and the large moment antiferromagnetism 
(LMAF). The first-order boundary separating these states is 
in accordance with the experimental datapii (open circles). 
For larger fields, only the HO survives, collapsing into a po- 
larized paramagnetic metal (PPM) b. Comparison of the the- 
oretical magnetic-field dependence of the normalized energy 
gap of the Q spin excitation in the HO, to the neutron scat- 
tering data^. In the low-field limit, the theoretical curve fits 
to the experimental one. The Q spin excitation is identified 
as an orbital antiferromagnetic resonance, c. Comparison of 
the temperature dependence of the normalized energy gap of 
the Qi spin excitation in the LMAF, to the normalized tem- 
perature dependence of the LMAF order parameter. The very 
good agreement implies a possible spin-exciton formation due 
to an anisotropic s-wave momentum structure of the LMAF 
order parameter, e. Predicted Polar Kerr effect due to the in- 
trinsic chirality of the proposed hidden order. The emergence 
of a chiral d-spin density wave below To — 17. 5K, implies a 
field-induced large non-zero Kerr angle, even above To- 

0.4,0,0) and the LMAF state Sq, of wave-vector Qq = 
(0,0,1) = Q = (1,0,0). Rotatmg from the Qi-dSDW 
state to the s-wave Q-LMAF, demands the presence of an 
idx^-y^ spin-0 orbital antiferromagnetic collective mode 
of momentum Qorh = (±0.4,0,0). We have constructed a 
Ginzburg-Landau theory for the vector order parameter 
n = {A2,S^), unifying HO and LMAPii. Although a 
phenomcnological model, it succeeds in explaining the 
p — B phase diagram (Figl3^) of the competing HO and 
LMAF, and the field and temperature dependence of the 
INS energy gaps. 

The HO exhibits an INS spectrum involving gapped 
spin excitations at Qi and Q. The first mode can 
be interpreted within the spin-cxciton scenario. The 
observed value of oJq^ — 4.5meV, is practically field 



independent^. According to the spin-exciton picture^^, 
the value uJq_^/2 = 2.25meV should be indicative^ of 
the energy gap of the HO. This value agrees well with 
A2 = 1.55meV, which is the dominant gap for the HO in 
our framework. On the other hand, the Q spin resonance 
originates due to the spin-0 orbital antiferromagnetic col- 
lective mode. At first sight, the singlet character of this 
mode and the momentum it carries do not seem to agree 
with the observed resonance. Nevertheless, the presence 
of the dSDW, provides the complementary symmetry 
breaking so that this collective mode indeed affects the 
INS response. FigI3]D presents a comparison of our theo- 
retical predictions compared to the experimental data^, 
focusing on the magnetic-field-dependence of the normal- 
ized energy gap a;^(2.2K,p = OGPa, i3^)/w^(2.2K,p = 
OGPa, OT). We observe that in the low- field regime, our 
results exactly fit the experimental measurements. 

In the LMAF state, a Bragg reflection appears at Q 
due to its pure AFM moment. Consequently, we expect 
suppresion of the Q resonant channel, which is indeed ex- 
perimentally verifiedii. However, the Qi resonance still 
persists in the LMAF state with a larger gajJ^. The 
existence of an orbital antiferromagnetic mode cannot 
provide this resonance, even if LMAF is present. A pos- 
sible explanation is that the LMAF order parameter has 
an unconventional s-wave momentum structure. In this 
case, the INS spectra would originate from the forma- 
tion of a spin-exciton at Qi, and the spin excitation 
gap would be proportional to the LMAF order parame- 
ter. Here we compare the theoretically predicted normal- 
ized order parameter temperature dependence, \n(T,p = 
0.67GPa,0T)|/|n(0K,p = 0.67GPa, 0T)|, to the nor- 
malized energy gap temperature dependence found in 
experiments^^. Indeed, we find a very good agreement 
(FiglS];) that supports the possibility of an unconven- 
tional s-wave LMAF order parameter. 

To detect the proposed chiral hidden order we pre- 
dict a Polar Kerr effect originating from the inherent 
time-reversal violation^!. This type of measurement is 
an ideal probe of the intrinsic chirality of the stated. 
More importantly, we expect that in our case the aris- 
ing Kerr angle , will demonstrate a magnetic field de- 
pendence when the sample is cooled down in the pres- 
ence of a static external magnetic field (FiglSli). We find 
that the Kerr angle probed by a Sagnac inteferometer 
(A = 1550nm), is 1?^ ~ 128Aimrad with Ai in mcV. 
The maximum value of A™°^ = 1.5meV, sets an upper 
bound -ff™"-^ ~ 190mrad. The magnitude of the Kerr 
angle demonstrates that the expected response is very 
large and easily detectable. The observance of a finite 
Polar Kerr signal could constitute a sharp signature of 
our chiral hidden order state, unveiling the mechanisms 
that lead to the intricate URu2Si2 phenomenology. 
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I. SUPPLEMENTARY MATERIAL 
II. MEAN-FIELD DECOUPLING AND ORDER PARAMETER SYMMETRIES 

Based on the phenomenology of URu2Si2, we have proposed that the chiral d-spin density wave state orders at 
an incommensurate wave vector Qi = (1 ± 0.4,0,0). Describing microscopically a chiral state with this ordering 
wave-vector, demands taking into account the main four Fermi lines that give rise to this nesting symmetry^. This 
is a very complex and unnecessary task, since our results do not depend on the microscopic details of the model. 
Our results concerning the anomalous Hall thermoelectric transport and the energy gap of the collective mode, are 
based solely on the topological content and the symmetry properties of this state. Moreover, the phase diagram 
and the metamagnetic transitions are also based on the general principles of itinerant metamagnetism, according 
to which this phenomenon originates from band crossing^. Under these conditions it is permissible to consider the 
existence of a nesting condition for the paramagnetic system e(fc + Qi) ~ ~e{k) = — [— 2t(cosA;a:+cosfcj,)] that is exact 
only for the commensurate ordering wave- vector (7r,7r). This condition, should stand for the real bands of URu2Si2, 
and within our approximate simple model we perform a simulation of this situation. Finally, we have considered 
Aj^(fe) = A]^ sinfca; sin fey, A2(fe) — A2(cosfc2; — cosky) with = ^1 1- 

To study the phase diagram of the chiral d-spin density wave state we consider that it is driven by an a inter-site 
Coulomb repulsion including nearest and next-nearest neighbours 

i^j i^j s,s' k,h' s,s' 

with s,s' =tii ^nd V{q) — 2V'{cosqx -f-cosQy) -I- 41/" cos cos . The momenta fc,fc' belong to the whole Brillouin 
zone {B.Z.) if Q is incommensurate and in the reduced Brillouin zone if Q is commensurate satisfying k -\- 2Q = k. 
The potential is separable providing V{k — k') = J2n^n fni^)fnif^') with n = 1,2,.., 8 corresponding to the form 
factors fn{k) : cosfc^^ ± cos sin ky , cos kx cos ky , sin kx sin ky , cos kx sin ky , sin kx cos ky with the driving 

potentials for n = 1, 2, 3, 4 and = AV" for the rest. Then within a mean-field treatment we introduce the 

staggered order parameters AQ^(fc), Aql'(fe), Aq^(A;), Aq (fc) as 



k' n 
fc' n 

and the corresponding irreducible fc-independent staggered order parameters 



\ss 

k 

The interaction decouples in the following way 



E {^Q ('^)4,tCfc+Q,T + ('^)44Cfc+Q4 + ^•'^■} + ^ E E - '^') (^Is^k+Q^s) (4'+Q,s'Cfc',, 

k k.k' s,s' 

E {a^q (fc)4,t^fc+Q,t + A^J-(fe)44Cfc+Q4 + /^--l + -E E (5) 



with V the volume of the corresponding B.Z,.. By introducing the spinor = (4 -|~ cJy. ^ 4+q t "^fc+Q i) ^^"-^ 
by using the isospin and spin Pauli matrices Tj,(Ti with i = 1,2,3 complemented by the related unit matri- 
ces Tq,(Tq, we furnish a representation of the 4x4 mean- field hamiltonian as Kronecker products of the form 
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Tfj, ® where [i^v = 0,1,2,3 (for simplicity wc omit ®). The most general chiral d-spin density wave state 
is defined as AQ(fc) = Aj^(fc) — iA2(fe). Moreover we introduce the kinetic term of the paramagnetic state 

(£(fe) - /i)ct ^Cfc., + (e(fc + Q) - Ai) along with a chemical potential ^. Putting together 

all these terms and using the fact that in the spin-triplet state Aq = Ag^ = — Aq, we obtain the quasiparticle 
Hamiltonian 



e{k)-fi Ai(fc) -iA2(A;) \ 

e{k)-fi -Ai(fc) + iA2(fc) 

A|(fc)+iA^(fc) e{k + Q)-fi 

\ -A|(fc) - iA5(fc) e(fc + Q)-/i / y 



k+Q 



Ai(fc) + A|(fc) Ai(fc)-At(fc). .A2(fc)-A;(fc) 

7, JTjCTa-i TiCTg 



A2(fc) + A^(fe) 



e(fe) + e(fe + Q) e(fc)-e(fc + Q) , ^ 



'2"3 



2 2 
and the elastic energy needed to build up the two density wave gaps 



(6) 



^A^=2«(^^ + -^j . (7) 

Wc now exploit the constraints that are imposed to the order parameters when the ordering vector Q is commensurate, 
satisfying fc + 2Q = k. The translation operation corresponds to the shifting k — > k + Q. Since k + 2Q = fe we 
have = 1 =^ ±1. This means that this discrete symmetry tranformation can divide a function into two irreducible 

representations, the periodic and the anti-periodic under Iq. This operation induces a transformation Iq in the spinor 
space, defined as 



iQ^^Q^^l'^O^fc- (8) 

An element in the spinor representation of the form '^\jl{k)'^ is invariant under translation, which is a symmetry 
transformation of the system, only if n(fc) transforms in the following way 



iQ^{k%=r^^^nk)T^^„. (9) 

In this manner we may label the elements of the spinor space with representations of the traslation operation. It 
is straightforward to show that the elements Tq ^a^ remain invariant under translation while the rest i.e. T2 3(7^ are 
anti-periodic changing their sign {v ~ 0, 1, 2, 3). This implies that for a chiral spin-triplet (and singlet) d-density wave 
we must have Aj(fc + Q) = Aj(fc) => A^ = A| and A2(A; + Q) = -/S.2{k) ^ IS.^ = -A2. Nevertheless, there is no 
such constraint when the ordering wave- vector is incommensurate like in our case. 



III. INTRINSIC ORBITAL MOMENT AND MAGNETIC FIELD COUPLING 

With these definitions we may write the quasiparticle hamiltonian in the condensed form 'Hg_p = 
J2k *fc {di^) ' ^ m} ^fc where we have introduced the vector g{k) = (^Ai{k)a^, A2{k)a.^, e(fc)crQ^ . Diagonalizing the 
above hamiltonian yields the eigenstates \^s,v{k)) with energy dispersions Es.v{k) ~ — /i -I- vE{k) and v = ±,s =t, 
E{k) = \g{k)\ = ^£'^{k) + |Aq(A;)P. The non trivial topological content of this Hamiltonian generates a non 
zero U{\) Berry connection^. A.s,y{k) = ($s_^(fc) | zV^ | ^s.iy{k)) and a corresponding Berry curvature defined as 
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^s,i^{k) = Vfc X A.sM{k) = iy{k)z, i.e. oriented along the z-axis. The detailed expression for the Berry curvature 
reads^"— 

^lAk) = '-TrJ^m . f W , . (10) 

^'"^ ' 2 ''\2E\kf^ ' \ dk, dky J j ^ ' 

with s, = ±. We observe that the Berry curvature is not depending on spin while it just changes sign when v = ±. 
The intrinsic orbital moment of this state is straighforward calculated using the definition^ 

m,,,(fc) = I X [Hik) - E^Jk)] I Vfc$,.,(fc)) = ^E{k)n,Jk) . (11) 

One may observe that the orbital moment is band independent. When we apply a magnetic field, the system interacts 
with both the Zeeman and the orbital moment mz{k) = eE{k)Vll _^{k)/h yielding the following interacting Hamilto- 

nian T-Lsik) = —{^B^Qa^ — mz{k)TQaQ)B. The inclusion of the magnetic coupling leads to the four field dependent 
eigenenergies Ef ,j{k) = — /i^s — rnz{k) — n + i'E{k) with z/, s = ± corresponding to the conduction and valence bands 
of the upper and lower spin sectors respectively. 

IV. FREE ENERGY AND SELF-CONSISTENCY EQUATIONS 

Having set up the microscopic description of the chiral d-spin density wave condensate, interacting with both Zeeman 
and the intrisic orbital moments with a perpendicular magnetic field, we may solve the self-consistency equations that 
derive from the minimization of the free-energy functional 



with t = 50meV, fj. = 0.69mcV, /i^ = 0.058mcV/T, a = 5A, V = 23.5meV and V" = 35.25mcV. Minimization of 
the free energy yields the following two self-consistency equations, that determine the order parameters of the chiral 
d-spin density wave^ 



A y^ ^ ..^ fUk) eaH^ ~ S{k) 



0^ A_^y _,A,iiA^_f^!i^S^A, 



aA, ' 2" fc^. I '"2-E{k) h 'E^{k) 



1 - 2 



\ E{k) ) 



[EfAk)] , (13) 



with ~ 4y", V2 = V , A-^ 2 ~ ^2 1' ^(k) ~ si^^ cos^ ky + sin^ ky and Up the Fermi-Dirac distribution. 

V. MAGNETIZATION AND THERMOELECTRIC RESPONSE 

As far as the magnetic and transport properties are concerned, we have used the expressions found in the theory 
of orbital magnetization^. Specifically, the magnetization is the summation of both spin and orbital contributions in 

.spin ' ZjOrb 

(14) 



the following way = -^Lpin + -^f 



^lor, - ^ E { (1 + ^"^.^(^)) ^^{k)np[El{k)] + '-^niAk)ksTln (l + e-^f..('^)/^-^) | . (15) 

The above expressions include the Berry phase correction to the electronic density of states, as this is apparent from 
the factor ( 1 + ^^^fif v{k) \ in the first equation. 
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As far as the transport properties are concerned, we have taken into account both topological [cr^y^top^ '^xy,top) ^^'^ 
quasiparticle contributions {(J^y^q-p, C(xy,q-p) calculating the Hall transport, while we have incorporated the Berry 
phase corrections in the scattering time and density states of the system. We have used Ts^v{k) = r = 3.3 • 10~^'^s. 
Based on experimental datai^ m* = 25me, and uJcTiB = IT) = 0.08, the estimated scattering time is of the order of 
10~^^s. In order to illustrate the chirality induced giant Ncrnst signal, we have on purpose selected a 3 times smaller 
relaxation time so to demonstrate that the possible enhancement of the quasiparticle contribution due to a large r, is 
not crucial for the anomalous Hall thermoelectricity. In finite magnetic fields, the initially constant relaxation time, 
becomes field and momentum dependent in the following manner 



The electric and thermoelectric coefHcients read 



B 

XX 



^xx 



^B 



xy,q-p 



.2^ i + ^mi^ik)\ n IK.WJ > (i7) 



<,top = E ^plEtimiAk) , (19) 



V 



- E U ] -AELm + In (l + e-^^C^)^-.-) i ntAk) , (20) 



k.s.u K \ rS 



= -- E i + p»g^(fc) [-IT [e.,^v->)—yy:^{k)^^^J=x,y, (21) 



2 



with units [a] = ea^/h and [a] = {kg/e)[a]. We have also introduced the field dependent velocities vfj^{k) — 
'V i^Ef ^{k). The Nernst signal and the resistivity are calculated from the following expressions 



Ey 

-dxT 



B B B B 

(J a — a (J 

xx^xy ^xx xy 
xx^ yy ' ^ xy xy 



rxx 



^ XX yy ' ^ xy^ xy 



(23) 



VI. GINZBURG-LANDAU THEORY OF THE MAGNETIC-FIELD-INDUCED CHIRAL D-SPIN 

DENSITY WAVE STATE 

We may construct a simple phenomenological model for the field-induced chiral hidden order proposed in this 
manuscript. The mumerical solution of the equations, points to the following Ginzburg-Landau picture for the chiral 
d-spin density wave order 

A? A2 A* 

F = a^^ + a{T-T,)^+P^-9^,^^B,. (24) 

The phenomenological constants a^^ji/^ ^^d the field B^^ are positive. The coefficient ensures that the d^y 
component is magnetic-field induced^. The important result discerned from the numerical analysis, concerns the sign 
of the orbital coupling g, that is also positive. With this convention, Aj^ 2 have the same sign. The mean-field solution 
for the Ai component yields 
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— = ^ a^A^ = . (25) 

One may 'integrate out' the induced Ai component and obtain an effective theory only in terms of the driving d^^-y^ 
order parameter A2. This is straightforward, yielding 



A2 a4 2 a2 a4 

Teff ^a{T-T,)^+/3^- ^AlBl = a[T-T, (SJ] -f + ■ (26) 

Notice, that since ai > 0, the quadratic magnetic-field-coupling leads to an enhanced critical temperature^ T^{B^) = 

+ ^B'^ for the dx^^y2 component and consequently for the chiral order. In the magnctic-field-induced chiral d-spin 
density wave, the A2 order parameter is driving the transition. 



VII. SO(2) MODEL FOR UNIFYING THE CHIRAL HIDDEN ORDER AND THE LARGE MOMENT 

ANTIFERROMAGNETIC STATE 



As we have previously shown, in the case of a magnetic-field-induced chiral d-spin density wave state, we may 
construct an effective theory solely in terms of the dominant d.j.2_y2 order parameter. Essentially, this is the order 
parameter of the chiral hidden order that we expect to demonstrate a degeneracy with the large moment antiferro- 
magnetic state. To handle the emergent SO (2) symmetry between the d-spin density wave component A2 of the chiral 
hidden order and the large moment antiferromagnetic state Sq, wc shall employ a simple S0(2) invariant Ginzburg- 

Landau model. The order parameters can be combined in a vector order parameter n ~ (A2,5q) as in reference^S. 
For a microscopic formulation based on operators, a more realistic model should be used. The phenomenology of 
URu2Si2 imposes the emergence of a first order transition from the hidden order to the antiferromagnetic state upon 
raising the pressure piiii^. At ambient pressure URu2Si2 5f-electrons condense into the HO state, implying that the 
latter should be promoted compared to antiferromagnetism. This calls for pressure dependent coefficients for the two 
order parameters that become equal at the critical pressure in which HO switches to LMAF. Since we study static 
order parameters our free energy functional is written 



T-T^^,ip,B,) 



^2 



r-T,_2(p,6J 



(5, 



1^2 



-13 



with 



T^,,{p,BJ = TM-gMBl. 



(27) 



(28) 
(29) 



We have introduced the positive coefficients a,/3 and the pressure dependent critical temperatures T^(j)),T^{p) and 
magnetic- field couplings go{p), g^{p)^ corresponding to the HO and the large moment antiferromagnetism, respectively. 
The critical temperatures and the coupling strengths will be considered to have the following pressure dependence 

Toip) = r„(o)+[T, -r„(o)] ^, (30) 

Pc 

^ V.T^-V^TM ^ TM^p , (31) 

Pi -Pc Pi- Pc 

5o(p) = v^ + UoP, (32) 

9aiP) = V^+^aP, (33) 



with Tj, the critical temperature at the critical pressure where the 50(2) symmetry is exact, T^{0) the critical 
temperature of the hidden order at zero pressure and T^{pi) the critical temperature of the antiferromagnetic phase 
for a pressure p^ > Pc- By setting n — (Aj, Sq), the free energy functional becomes 
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1=1,2 ^ 

Minimizing the free energy functional, provides the equations 



2 



i2{|np-f [r,,2(p,SJ-T]}=0. 



(34) 

(35) 
(36) 



The above system of equations, demonstrates that the hidden order and the antifcrromagnetic states are competing 
and are separated by first-order Une phase boundaries. 



VIII. INELASTIC NEUTRON SCATTERING ENERGY GAP-MASS OF THE SPIN-0 ORBITAL 

ANTIFERROMAGNETIC COLLECTIVE MODE 

The presence of the orbital antifcrromagnetic mode with Qorb = (iO.4, 0, 0) momentum, may be recorded indirectly 
through the pressure and field dependence of the energy gap of the spin excitations. When the SO (2) symmetry is 
broken, this mode acquires a massiSiHiil that depends on pressure and field. This mass corresponds to the energy gap 
of the spin excitations. When p ^ pc and = 0, the SO (2) symmetry emerges and the mode is masslcss. We may 
calculate the arising mass of this mode in terms of the phenomenological model of the previous section, developed for 
the order parameter n. We perform a % rotation — >■ rii cosx + '^2 sinx and rtj cosx — sinx, that yields 

the X dependent free energy functional 



nl cos^ X + '^2 si^^^ X + '^1^2 2x 



.2 ,-„2 



T-T^Jp,B,) 



cos'^ X + ^1 sin^ X ~ '^1^2 ^™ 2x 



n 



(37) 



By expanding J-{x) in powers of x about its zero mean-field value, we obtain an effective free energy functional 
descriding the static part of the orbital antifcrromagnetic fluctuations. The linear term in x, necessary disappears 



dx 



dT _ a 

^ ~ 2 
a 

+ 2 

— — ( 



T - T^ iip, B^) y-n\ sin 2x + nl sin 2x + 2n^n2 cos 2x 
T -T^ 2{p, B^) (-nl sin2x + nf sin2x - 271^712 cos2x 



0. 



(38) 
(39) 



since on one hand 77^^ and 772 never coexist for p ^ p^, while on the other, their temperatures become equal exactly at 
Pj,. The second term in the expansion is associated with the mass of the collective mode. 



dx^ 



a 



T - T^ -^{p, B^) (^-n\ cos2x -t- 773 cos2x - 277^772 sin2x 
T — T^ 2{Pi ^z) ( "^^2 cos2x -I- 71^ cos2x -I- 277^^772 sin2x 



9oip)Bl 



(40) 
(41) 



a{n\ - 772) 



PlTc - PcTa iPl) , T^ iPi)~T^ 



Pi ~Pc 



Pi "Pc 



■P " 9a{p)^l 



(42) 
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The mass of the collective mode m is given by the following definition Fmass = \m?x^- By taking the second 
derivative, we finally obtain 



2 / 2 

m = a [rii 



Pl -Pc 



[Oaip) - 9o{p)] Bl 



(43) 



Once again, since the order parameters ni,n2 never coexist, we have — n| = |n|^ = in the HO and n1 — 



-712 ill the LMAF. We must remark that the mass of the mode m should be related with the resonant 



excitation at q = Q in the HO. 



Pl -Pc 



Pc/ 



(44) 



The contribution of the kinetic energy of the mode to the free energy functional, is expressed as —^pu'^x^- P is a 
coefficient that would arise if we had calculated the fluctuating part of the collective mode. Such a treatment would 
demand a microscopic model that could properly take into account the main Fermi surface sheets of this compound. 
This coefficient is defined in meV and shall be considered constant. The fact that we do not exactly know this 
parameter, prevents us from obtaining the exact energy gap. However, we may infer the normalized field and tem- 
perature dependence of the energy gap of a spin excitation carrying momentum q, defined as Wq(T,p, B^)/u!q{T,p, 0) 
and uj'^{T,p,B^)/uj'^{0,p,B^), respectively. 



^"g{T,p,BJ 

^°{T,p,0) 



n{T,p,B, 



niT,p,0) 



\ 



Pc{pi-Pc^ 9a{p)-go{p) 



(pc - p) {pi [To (0) - Te] - \To (0) - Ta (Pi) 



Bi 



1 P^Pc 



(45) 



In order to continue, we use a number of experimental data so to identify the HO and LMAF critical temperatures 
and magnetic couplings. The specific pressure dependence that we have considered for the critical temperatures, will 
give a phase boundary that has a vertical line in the p — B phase diagram. However, experimental evidenceii show 
that the first-order phase boundary is a straight line only for temperatures below ^8K. In this case p^ = 0.5GPa, ~ 
17. 8K—. This choice of p^, restricts the validity of our phenomenological theory to low temperatures. However, the 
neutron scattering measurements are usually performed below 8K. Consequently, our model works well in this region 
and we may compare our calculations to the experimental data. Specifically we use: (OGPa) = 17. 5K at Bz = ^^'^^ , 
= 2GPa, T^iPi) = 19. 5K at B^ = Oii, *'°(OGPa) ^ 40T at T = OK-i^ and that *'°(0.72GPa) = 15T at 
T ~ 16K— . Moreover, we have used the experimentally extracted dependence ujq{2K, OGPa, Bz)/ujq{2K, OGPa, OT) = 

\/l + 0.005 • B'^, for the Q = (1,0,0) resonance in the HO. Finally, from Ref.— , we have used that for pressure 
p = l.lGPa and temperature T=OK, the LMAF state undergoes a first-order transition to the HO state for Bz = 12T. 
At this critical point of the phase diagram, the orbital antiferromagnetic collective mode becomes massless, providing 
zero energy gap to the spin excitation spectrum. The arising condition from this experimental feedback may interpreted 
as a;q(OK, l.lGPa, 12T) = 0. The aforementioned conditions define in a unique manner the parameters of our theory 
providing the dependences: T^{p) = 17.5 + 0.6p(K), T^{p) = 17.2 + 1.13p(K), g^{p) = 0.011 - 0.003p(K/T2) and 
9a(.P) = 0-012 — 0.019p(K/T^). The normalized magnetic-field-dependence of the energy spectrum reads 



^°(T,p,gJ _ n(T,p,g, 
o^°{T,p,0) 

Finally, for the two resonances we have: 



n{T,p,0) 



Up 



0.5 



(O-OlBzY 



P'^Pc 



(46) 



^"QiT,P,Bz) 
^%{T.P.Q) 

^°qAT,p,Bz) 
^'qST,P,0) 



^ 0.011 - 0.003p 
17.5 + 0.6p-r ^ 



0.012 - 0.0019p 

1 — B 

17.2+1.13p-r ' 



1 + ?^(0.0W 




0.5 -p 



25 -I- 12p 
0.5 -p 



(O.OlBzy 



T < 8K (HO) , (47) 
T < 8K (LMAF) . (48) 
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FIG. 4: Interpolation of experimental data taken froii>ii, concerning the temperature dependence of the energy gap of the Qi 
spin excitation, observed in the large antiferromagnetic state. 



We have obtained the p — B phase diagram at T=OK, by retrieving all the combinations of pressure and field that 
reduce the above expressions to zero. These values have provided the first-order boundary of the HO and the LMAF 
and the transition from the HO to a polarized paramagnetic metal (PPM). In addition, we used Eq|47l to compare 
the theoretical outcome for cjq(T, OGPa, S^)/wq(T, OGPa, 0) with the corresponding experimental valuesii. 

The normalized temperature dependence of the spin excitation energy gap becomes 



1 - 



T 



17.5 + Q.6p - (0.011 - 0.003p)i32 



1 - 



17.2 + 1.13p - (0.012 - 0.0019p)i32 



T <8K (HO) , 
T < 8K (LMAF) 



(49) 
(50) 



We oberve that in this case all the parameters are fixed. The conditions that we have used to define the values of 
these parameters were uncorrelated to any information concerning the temperature dependence of the energy gap of 
any of the possible resonances. Any theoretical calculation concerning the temperature dependence will constitute a 
prediction of the S0(2) theory. 

We go on with proposing an explanation for the observed Qi resonance in the LMAF state. The spin-0 orbital 
antiferromagnetic mode cannot couple to a Qi modulated external field, B^^, even in the presence of the LMAF. 
The reason is its dg.2_y2 momentum structure. However, a resonant excitation is possible to arise in the case of an 
anisotropic s-wave order parameter for the LMAF, that could provide appropriate coherence factors for the formation 
of a spin-exciton at Qi. In order to test this possibility, we have made a comparison to experimental measurements. 
For this task, we have used the data presented in FigHl that were taken from Refiii and concern the temperature 
dependence of the energy gap of this mode. Since, there were not many available data below T=8K, which is the 
limit of the validity of our theory, we have interpolated the experimental data. We have compared the interpolated 
data uJq{T)/u!q{0) ~ (8.28 — . 28r) /8 . 28 to the normalized energy gap of EqlSD) The theoretical and the experimental 
curves are in very good agreement indicating that the normalized energy gap of the resonance mode scales with 
|n(T,p, 0)|/|n(0,p, 0)|. This clearly implies that the spin excitation develops a gap proportional to the LMAF order 
parameter, supporting our proposal. 



IX. POLAR KERR EFFECT 

The Kerr angloi^ri^i is defined as 

K ~ I 9 TTTT ' ' 

where n is the refractive index, d is the c-axis thickness of the sample, a = 5A the in plane lattice constant and 
a'^y the imaginary part of the dynamic Hall conductivity. In the low field regime {B < 20T), the dynamic Hall 
conductivity is dominated by the topological contribution. For a Sagnac inteferometer— huj = O.SeV = 6420cm~^ 
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and A = 1550nm, while the plasma frequency of URu2Si2 is hiOp = 2131cm ^ ~ 266meV— . The refractive index in 

\2 



this case is defined as n = y 1 — {ujp/uj) . Since uj » Up we may use the approximation ti ~ 1 and n — 1 = — {ujp/uj) 

For the calculations of the Polar Kerr effect in the presence of small external static magnetic fields, only the topological 
dynamic Hall conductivity is important in our case, defined as 

—iuj 

with n^^(<7,icj^) the current-current charge polarization function originating from the bubble diagram 

2 

n-(za;„ g) ^ j {^(*^- (^^ ^ + q)Tp g{ik^ + iLU^,k + q) r^-(fc + q, fc)r,} (53) 

and g{ik„, k) is the bare matrix Matsubara Green's function 

1 — iJ.g(J^ — m^(k) —1.1 + g{k) ■ T 



g{ik^,k) 



ik„ - /.(^(Tg - mz{k) - /i - g{k) ■ t (ifc„ - ^^cTg - mz(fc) - fi)^ - E'^{k) 
1 f E{k)+g{k)-T E{k)-g{k)-T 



2E{k) 1 j/jj^ _ Mfifa ~ nizik) — /J — £'(''') i^n — Ms'^'s ^ niz{k) — /J + £-(&) 



(54) 



and the vertex functions r^{k,k + q) are defined as the coefficients of the fermion-gauge field coupling written 

in the form J2q — J2q k^k+q^i(^ + 9; ^)'r/i^fc^*(~'Z)- Consequently they are known if we find the 

fermionic currents. We start by the equation of continuity of the electric charge Qez = / drp{r) in momentum space 

p{g) ~ig' J'^ig) — O, where we have introduced the electric charge density p(q) = — eX^fc^I+q^fc- The corresponding 
electric current can be derived from the equation of continuity as the limit 



r[q) = -I lim V, p{q) = lim [ 7^, p(q) ] = -e J] vj/t (fc)*^ (55) 

''^ k 

with V{k) — V^'H(fc) the velocity defined in spinor space. Consequently r°(fc + q, fc) = and T.^ {k + q,k) — '^^^''^ . 
Straightforward calculation of the polarization tensor yields the Hall topological conductivity 



e2 ^ 4i?2(fe) Up [Ef,,{k)] nijk) 



hv [hLO + 27/ - 2E{k)] [huj + i-q + 2E{k)\ ' 
For w > the Hall conductivity has a contribution from only the first pole and one may write 



e2 ^ AE^{k)n^\Ef(k)\D.l(k) ( V ^ 
ktu /iw + M/ + 2£;(fe) \huj-2E{k) ^ ^ '^j 

where V denotes the principal value of the sum. The imaginary part of the conductivity is given as 



(56) 



where due to the delta function we have the simplification Ef jj{k) = —jj— [i-LBS + mz{k)]B + i'fkij/2. Typically, haj is of 
the order of eV, definitely larger than the energy scales considered here. This means that np ^.(fc)] ~ np[i^huj/2]. 
For low fields, the physics may be well described by considering an expansion around the feg — (fj §) point. In this 
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case, ni^{k) ~ 2iyaHAiA2/E^{k) = Wi^aHAi A2 / {hiu)^ , E{k) = ^ Al + Al{5k^ - Skyf + {2t)^{5k^ + 5kyY, with 
Sk = k-ko and i EkeB.z. = dE E/[2^(4tA2)]. 



cr^ (w) = — r 7 a ~ -ri— a (57) 



The Kerr angle is equal to 



~ Ijaa^ u n[n'^ — 1) he a noj 



nc a hUp 

By taking into account the value of the fine structure constant /he = 1/137 and by setting d ~ 2nni, we obtain 
a rough estimation of the Kerr angle as -d^^ ~ 128Ainirad with Ai in meV. The maximum value of Ai is equal to 
1.5meV, which defines an upper bound i?^ < 190mrad. 
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